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SECTION  1 


INTRODUCTION 


1.1  SUMMARY 

This  report  summarizes  the  results  of  an  investigation  of  bi-directionally 
reinforced  :oinpoalte  flywheels.  The  investigation  focuses  on  the  hoop-radial 
type  of  bi-directional  configuration  developed  by  Avco  and  includes  both 
analytical  and  experimental  work.  This  work  is  divided  into  four  major 
portions : 

^ Analytical  Assessment  of  the  Effect  of  Degraded  Material  Properties 
(Section  2) 

• Sensitivity  Studies  (Section  3} 

• Analysis  of  Static  Simulation  Test  (Section  4) 

• Testing  for  Degraded  Property  Effects  (Section  5) 

Some  conclusions  from  this  work  are:  (1)  decreases  of  in-plane  shear 
modulus  of  up  to  90  percent  have  small  effect  on  the  maximum  stresses, 

(2)  changes  in  hoop  thickness  or  fiber  volume  ratio  can  cause  changes  of 
similar  magnitude  in  theoretical  flywheel  efficiency,  (3)  static  simulation 
tests  are  feasible  but  must  Include  substantial  numbers  of  radlals  to  prevent 
premature  shear  failures,  and  (4)  degraded  property  effects  can  be  measured 
experimentally  and  substantial  reductions  in  moduli  and  shear  stress  allow- 
ables are  found.  These  and  other  results  are  discussed  in  greater  detail  in 
Section  6. 

1.2  BACKGROUND  AND  PRELIMINARY  DEVELOPMENT 

Bi-directionally  reinforced  flywheels  offer  potential  advantages  over 
the  more  conventional  composite  ring  designs,  such  as  greater  energy  per  unit 
volume  and  a larger  number  of  geometric  parameters  to  apply  to  the  design  of 
hub  connections.  In  addition,  theoretical  analyses  show  that  the  stored 
kinetic  energy  per  unit  weight  of  a bi-directionally  reinforced  wheel  can  be 
within  a few  percent  of  that  obtained  through  ring  designs. 

The  kinetic  energy  per  unit  weight  of  a composite  flywheel  composed  of 
one  material  can  be  written  in  the  following  form. 

Kinetic  Energy  ^ a 
Weight  ® pg 


where 


Kg  * shape  factor 

a « allowable  tensile  stress  of  a 1-D  ply  in  the  0°  direction 

pg  weight /unit  volume  of  the  wheel. 

For  a given  material  a and  pg  are  known  and  the  primary  problem  is  to 
configure  the  wheel  to  achieve  the  highest  possible  value  of  Kg,  Analyses* 
based  on  the  neglecting  of  the  very  low  90°  tensile  strengths  indicate 
that  the  maximum  Kg  obtainable  for  a composite  is  .5.  (For  a simple 
ring  Kg  ■ .5.)  Design  experience  at  Avco  indicates  that  some  deviation  from 
the  optimum  is  necessary  to  prevent  excessive  hoop-radial  interlaminar  shear 
stresses.  As  a result  the  shape  factor  Kg  must  be  slightly  reduced  from  .5 
to  .474.  It  should  be  noted  that  since  the  formula  for  kinetic  energy  per 
unit  weight  involves  only  the  strength  in  the  0°  direction,  the  formula  ap- 
plies to  a flywheel  which  has  been  designed  to  avoid  premature  failures  as- 
sociated with  other  stress  components. 

The  Avco  bi-directional  flywheel  design  consists  of  a stacked  sequence 
of  hoop  and  radial  layers.  The  hoop  layers  are  approximately  .017-lnch  thick 
and  essentially  axlsymmetrlc.  The  radial  layers  each  contain  144  radial 
members  .008-lnch  thick  and  of  tapered  width  to  provide  optimum  efficiency. 
The  material  is  Kevlar  49.  A simplified  schematic  drawing  of  the  Avco  design 
is  shown  in  Figure  1-1.  This  schematic  drawing  shows  only  16  of  the  144 
radial  members  present  in  the  actual  design.  This  design**  is  used  as  a 
reference  point  for  the  investigations  carried  out  for  this  report.  Conse- 
quently, it  is  appropriate  to  outline  and  define  the  basic  parameters  that 
controlled  the  reference  design. 

At  the  onset  we  Introduce  appropriate  nondimens ional  stress  nomenclature 
as  follows : 

. a 

pfl2  r2 

where 

a m stress 

a > nondimens ional  stress 

p ■ mass  density  of  a ply 


•Johnson.  D.  E..  Maximum  Shape  Factor  for  Cor..po5ltB  Flywheel  Avco  Technical  Release  K400'T-322-A,  August  27, 1975. 
••Sapowith,  A.,  Phase  I Interim  Technical  Report.  Avco  AVSD-0189-76-CR,  June  18, 1976. 
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Figur*  1-1.  Schematic  of  typical  bi-directionally  reinforced  flywheel. 
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fl  » angular  velocity  in  radians/sec. 

R « outside  radius  of  the  flywheel 

The  nondimens ionalization  symbol  ' is  applied  to  various  stress  components. 

The  nondimensional  aspect  of  this  problem  implies  that  the  results  from  these 
analyses  can  be  applied  to  a whole  family  of  wheels  of  varying  outside 
radius  — provided  the  thicknesses  of  the  individual  plys  are  scaled  propor- 
tionately with  the  outside  radius. 

The  dominant  feature  of  an  efficient  bi-directional  flywheel  is  a biaxial 
state  of  stress  with  approximately  equal  stresses  in  the  hoop  and  radial 
directions.  For  Kevlar  the  0°  and  90°  failure  strains,  e and  (gn  fonn  the 
following  ratio  ° 


'90  0.22% 


This  implies  that  the  transverse  (90°)  properties  of  the  Kevlar  will  be  de- 
graded, and  crazing  will  occur  well  before  the  0°  failure  strains  are  reached. 
Because  of  this  reason  the  radial  configuration  for  the  reference  Kevlar  fly- 
wheel was  designed  on  the  basis  of  Eqg  ■ 0 in  each  hoop  and  radial  ply.  To 
implement  the  iterations  required  during  the  design  process  a "smeared  out" 
axlsymmetric  finite  element  model  was  used  to  obtain  the  optimally  shaped 
radial  reinforcements  that  yielded  the  highest  shape  factor  K3.  The  non- 
dimensional  stresses  derived  from  the  axlsymmetric  model  are  shown  in 
Figure  1-2.  This  figure  clearly  shows  the  dominant  features  of  the  problem: 
nearly  biaxial  ( cr,  • oq)  stresses  over  the  interior  of  the  wheel  and  very 
much  smaller  shear  stresses  and  (see  separate  scale  on  right  side  of 
the  figure).  Unfortunately  the  shear  allowables  are  correspondingly  smaller, 
as  may  be  seen  from  Table  1-1  which  summarizes  the  maxima  of  the  stresses 
shown  in  Figure  1-2  together  with  the  corresponding  stresses  from  a more 
complex  two-dimensional  finite  element  model  of  the  same  problem.  The  2-D 
model  is  thus  used  as  a check  on  the  simpler  equivalent  axlsymmetric  model 
which  is  so  well  suited  for  the  parametric  variations  of  radial  geometry  re- 
quired during  design. 


TABLE  1-1.  STRESSES  IN  REFERENCE  FLYWHEEL  DESIGN 
FOR  Ego  * 0 (DEGRADED) 


Equivalent  axlsymmetric 

Type  of  analysis 

Stress 

2-D  Finite  element 

“ Allowables* 

At  fi  • 5685  At  n - 5630 

Nondimensional  radians/sec  Nondimensional  radians/sec 


"r 

.516 

200  ksi 

.514 

200  ksi 

188-211  ksi 

.516 

200  ksi 

.530 

200  ksi 

188-211  ksi 

.00401 

1.55  ksi 

.00355 

1.34  ksi 

1.45-7.9  ksi 

.00982 

3.81  ksi 

.01032 

3.90  ksi 

7.9  ksi 

*Based  on  Avco  tests  and  duPont  report  "Characteristics  and  Uses  of  Kevlar  49 
High  Modulus  Organic  Fiber,"  by  R.  L.  Hunter,  revised  9/18/73. 
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Material  Properties 


Unless  otherwise  specified,  the  following  Kevlar  49  mechanical  properties 
are  used  in  this  report: 


Property 

Equivalent 

nomenclature 

Values 

Undegraded 

used 

Degraded 

El 

Eq  - E 

12  X 10^  psi 

12  X 10 6 psi 

E2  • E3 

Ego 

0.8  X 106 

0 

PI2  - 1^13 

.34 

0 

‘'21  ■ ‘'31 

.0227 

0 

‘'32 

.48 

0 

Gi2 

G 

0.3  X 106  psi 

. 3 X 10^  psi 
or  reduced* 

These  properties  are  derived  from  the  properties  given  by  duPont,**  As 
pointed  out  in  the  previous  section  the  degraded  properties  were  used  in  the 
flywheel  design  because  of  the  very  low  transverse  (90°)  strain  capability  of 
Kevlar.  It  should  also  be  noted  that  the  linear  elastic  nature  of  the  prob- 
lem is  such  that  the  resulting  stresses  are  unchanged  if  the  moduli  ^2* 

E3  and  G are  multiplied  by  the  same  scalar  constant. 

Comparison  of  Finite  Element  Models 

Three  different  types  of  finite  element  models  were  used  in  this  report: 

1.  1-D  Axlsymmetric  Model 

2.  2-D  Model  including  Discrete  Radials 

3.  2-D  Model  using  "Smeared  Out"  Properties 

Each  of  these  three  models  assumes  plane  stress.  Solutions  from  the  1-D 
axlsymmetric  model  are  functions  of  the  radial  coordinate,  r,  only.  Solu- 
tions of  the  2-D  models  are  functions  of  both  r and  6 . 

Examples  of  the  finite  element  grids  used  in  these  three  models  are  shown 
in  Figures  1-3  through  1-5.  In  the  2-D  model  with  discrete  radials  shown  in 
Figure  1-4,  the  nodes  in  the  hoop-radial  overlap  region  each  have  one  set  of 
planar  displacements,  l.e.,  in  spite  of  the  overlapping  the  model  remains 


•All  investigation  in  which  G112  is  reduced  below  this  value  is  presented  in  Section  2 of  this  report. 

••Hunter,  R.  L , "Characteristics  and  Uses  of  Kevlar  49  High  Modulus  Organic  Fiber,"  duPont  Report,  revised  9/18/73. 
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Figure  1-6.  Finite  element  grid  used  for  asymmetric  property  variations. 
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planar.  This  model  uses  2488  nodes  and  3410  elements  to  represent  the  very 
slender  angular  section  shown  in  Figure  1-4. 

Most  of  the  work  was  done  using  the  first  two  models,  the  last  being 
mainly  used  to  determine  the  effects  of  gross  asymmetries  and  imbalance. 

Each  of  the  three  finite  element  models  produces  the  stresses  in  the 
radials  and  the  a0$  stresses  in  the  hoops  directly.  In  the  2-D  models  the 
in-plane  shear  stresses  are  also  produced  directly.  Because  of  the  axi- 
symmetry  of  the  1-D  finite  element  model,  it  gives  o^q  ■ * 0*  Neverthe- 

less, it  is  possible  to  start  with  the  oqq  from  the  1-D  modal  and  compute, 
via  static  equilibrium,  the  stresses  that  would  occur  if  one  had  a set 
of  discrete  radial  reinforcements  as  shown  in  Figure  1-1.  The  essence  of 
this  calculation  is  shown  in  Figure  4-5  and  is  based  on  the  previously  cited 
damaged  properties.  None  of  the  three  finite  element  models  used  in  this 
work  gives  the  interlaminar  stress  directly.  Consequently,  the 
stresses  were  obtained  from  the  other  stresses  by  using  static  equilibrium. 

The  method  applied  to  the  1-D  model  is  very  similar  to  that  used  for  obtaining 
OfQ  and  is  outlined  in  Figure  4-5.  The  stresses  for  the  2-D  model  were 
also  derived  using  static  equilibrium:  In  this  case  the  net  force  on  a small 
element  of  each  layer  (hoop  and  radial)  were  computed  and  used  to  determine 
the  force  transferred  between  the  layers  bv  the  interlaminar  shear  stress. 

Because  the  state  of  stress  in  the  bi-directionally  reinforced  wheel  is 
almost  axisymmetric,  the  stress  distributions  presented  in  this  report  are 
generally  plotted  as  functions  of  the  radial  coordinate  r.  The  relatively 
small  deviations  from  axlsyrametry,  when  they  appear,  are  indicated  by  showing 
the  maximum  and  minimum  values  with  respect  to  changes  in  the  6 (circumfer- 
ential) coordinate.  In  figures  such  as  1-8  or  2-3,  the  deviation  from  axi- 
symmetry  is  sufficiently  great  so  that  two  separate  dashed  curves  are  used 
to  show  the  maximum  and  minimum  values. 

Because  these  several  different  models  were  used  to  carry  out  the  work 
presented  in  this  report,  it  is  appropriate  to  pause  and  obtain  a check  of 
the  accuracy  by  comparing  results  obtained  by  the  different  models  for  the 
original  flywheel  design.  This  was  carried  out  for  two  different  cases,  one 
using  the  undegraded  properties  shown  in  the  table  and  one  using  the  degraded 
properties  (with  G ■ .3  x 10^). 

Figures  1-6  through  1-9  show  that  for  the  hoop  and  radial  stresses  the 
two  models  give  almost  the  same  results  for  each  casj.  the  difference  between 
them  being  of  the  same  order  as  the  circumferential  variation  (see  oqq  ) within 
the  2-D  model  itself.  Similar  agreement  is  seen  in  Figure  1-10  for  the 
shear  stresses.  The  stresses  shown  in  Figure  1-11,  however,  show  a con- 
siderable difference  between  the  average  values  derived  from  the  axi- 
symmetric  model  and  those  from  the  2-D  model.  This  latter  difference  is  not 
accountable  from  the  circumferential  variations  in  the  2-D  model  but  appears 
to  be  caused  by  a fundamental  difference  between  the  two  models.  Because  of 
the  planar  nature  of  the  2-D  model  it  gives  a that  does  not  approach  zero 
at  the  edge  of  a radial.  This  can  violate  the  shear-free  condition  actually 
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Iwo  dimensional  model 
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Coarse  grid  two-dimensional 
Model  shown  in  Figure  1-B 
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Figure  1-9.  Original  flywheel  with  undegraded  properties,  comparison  of  radial  stresses  in 


fiber  direction  along  centerline  of  radial  reinforcement.* 
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Two-dimensional  model 
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SECTION  2 


ANALYTICAL  ASSESSMENT  OF  EFFECT  OF  DEGRADED  PR,OPERTIES 


2.1  INTRODUCTION 

In  this  section  we  investigate  the  effect  of  degraded  material  properties 
caused  by  the  high  biaxial  stresses  in  the  flywheel.  The  analytical  results 
are  obtained  by  using  the  two-dimensional  finite  element  grid  shown  in  Figure 
1-4  and  described  in  Section  1, 

2.2  PRELIMINARY  INVESTIGATION;  UNDEGRADED  VERSUS  DEGRADED 

As  previously  mentioned,  the  flywheel  is  designed  using  the  E90  ■ 0, 

■ 0 degraded  properties  because  of  the  low  transverse  (90°)  strain  allowable 
of  Kevlar.  The  resulting  flywheel  has  thus  been  optimized  on  the  basis  of  the 
degraded  properties.  If  the  same  configuration  is  analyzed  using  undegraded 
properties,  the  stress  distributions  are  as  shown  in  Figures  2-1  and  2-2.  It 
is  immediately  apparent  that  the  design  was  optimized  for  the  degraded  con- 
dition because  the  maximum  stresses  are  substantially  larger  in  the  undegraded 
design,  as  summarized  in  the  following  table; 


RESULTS  FOR  2-0  FINITE  ELEMENT  ANALYSIS 


Of 

max 

Undegraded 

Degraded 

max 

r2 

.514 

.605 

X.. 

pn2  r2 

.530 

.470 

Kg  (2-D  Analysis) 

.461 

.404 

This  table  also  shows  the  shape  factor  computed  for  these  two  cases  and 
indicates  a 12  percent  lower  Kg  (and  hence  lower  energy  density)  for  the  un- 
degraded condition  than  for  the  degraded  condition  used  to  optimize  the 
design. 

2.3  EFFECT  OF  DEGRADED  SHEAR  MODULUS 

A series  of  two-dimensional  finite  element  runs  were  made  to  determine 
the  effect  of  reducing  the  shear  modulus  G for  the  degraded  properties  shown 
in  Section  1.  The  results  shown  in  Table  2-1,  are  somewhat  startling.  They 


Indicate  that  the  maximum  stresses  and  shape  factor  are  almost  unaffected 
by  reductions  In  the  shear  modulus  dovra  to  one-tenth  of  the  orglnal.  The 
distributions  of  the  stresses  shown  In  Table  2-1  are  given  In  Figures  2-3 
through  2-5  for  the  two  cases  G/GQj.j^g  « 1 and  G/Gorig  “ .!•  Although  the 
figures  show  very  little  change  In  the  maximum  stresses  with  the  reduction  of 
G/Ggrlg  to  .1,  they  do  Indicate  the  beginnings  of  changes  In  the  stress  pat- 
tern. Figure  2-3  shows  that  as  the  shear  modulus  Is  lowered,  the  distribution 
of  stresses  In  the  circumferential  direction  across  a radial  becomes  In- 
creasingly nonunlfom.  At  G/Ggrlg  “ •!»  as  shown  In  Table  2-1,  the  nonunl- 
formlty  Is  already  sufficiently  great  for  the  radial  stresses  to  become  the 
maximum  stresses  In  the  wheel.  Figure  2-A  shows  a similar  circumferential 
nonuniformity  In  the  hoop  stresses  near  the  outside  of  the  wheel  as  the  shear 
modulus  Is  reduced.  These  trends  would  Indicate  that  the  and  stresses 
will  Increase  as  G Is  reduced  below  .1. 

TABLE  2-1.  EFFECT  OF  DEGRADED  SHEAR  MODULUS 
(FROM  2-D  FINITE  ELEMENT  ANALYSIS) 


G/^orlginal 

1.0 

.5 

.1 

\ax 

pn2  r2 

.514 

.517 

.526 

X.. 

p r2 

.530 

.529 

.523 

.00355 

.00356 

.00354 

^tO 

pil^  r2 

.01032 

.00958 

.00947 

Ks 

.461 

.462 

.465 

Implications  of  this  result  are  as  follows.  First,  within  the  range  of 
G degradations  considered  it  indicates  that  the  dominant  effect  on  flywheel 
efficiency  will  come  from  reduced  shear  stress  allowables  rather  than  from 
the  reduced  modulus.  The  differing  effects  of  shear  modulus  and  shear  allowa- 
bles are  Illustrated,  by  example,  in  Section  2.4.  Second,  the  results  Indicate 
that  G/Gorlg  niust  be  reduced  below  .1  to  alter  the  stresses  a-  and  suffi- 
ciently to  effect  the  overall  flywheel  shape  factor  Kg.  Further  evaluation  of 
extremely  degraded  G cases  Is  being  continued  at  the  present  time. 
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2.4  CORRELATION  OF  PREDICTIONS  WITH  MEASURED  DEGRADATION 

We  now  correlate  the  analytical  predictions  with  the  property  degradation 
measured  experimentally  and  reported  in  Section  5.  The  experimentally  deter- 
mined shear  modulus  ratio  G/Goriginal  is  *^25  and  falls  within  the  range 
covered  by  Table  2-1.  This  indicates  that  the  wheel  stresses  are  essentially 
unaffected.  The  reduction  in  shear  strength  to  .728  of  the  undegraded  strength, 
however,  indicates  lower  margins  of  safety  with  respect  to  shear  failure  for 
the  degraded  case.  For  the  undegraded  case  the  margin  of  safety  on  the 
stress,  computed  from  the  values  shown  in  Table  1-1,  is  (7. 9/3. 9) -1  = 1.02. 

For  the  experimentally  measured  degraded  shear  strength  it  is  (5.75/3.9)-!  * 

.47,  so  that  the  margin  is  reduced  roughly  by  half  but  nevertheless  remains 
positive. 

The  greater  sensitivity  of  the  margins  to  shear  strength  rather  than 
shear  modulus  can  be  further  illustrated  by  considering  two  hypothetical  cases 
and  the  resulting  margins  of  safety  as  follows: 


Case 

For  G 

Percent  of  undegraded  values 

For  <>12 
allowable 

Margin  of 
safety 

Case  1 

m% 

40% 

-.19 

Case  2 

40% 

100% 

+1.0 

Undegraded  Case 

100% 

100% 

+1.0 

The  table  Indicates  greater  sensitivity  of  margins  to  the  <^12  shear 
allowable.  However,  it  should  be  pointed  out  that  this  computation  of  mar- 
gins based  on  Otd  and  the  shear  allowable  may  be  conservative  because  even 
if  the  0^12  shear  strength  is  lost,  the  performance  of  the  wheel  may  not  be 
degraded  due  to  catenary-like  action  of  the  degraded  hoop  reinforcements. 

Because  the  experimental  results  apply  to  the  case  of  52  percent  of  the 
ultimate  0°  stress,  it  is  possible  that  the  degradation  of  In-plane  shear 
strength  may  limit  the  flywheel  performance  at  higher  percentages  of  the 
ultimate  load.  Such  a problem,  if  it  occurs,  could  be  ameliorated  by  in- 
creasing the  number  of  radlals  and  thereby  reducing  the  stresses. 


-34- 


SECTION  3 


SENSITIVITY  STUDIES 


3.1  INTRODUCTION 

Sensitivity  studies  were  carried  out  to  evaluate  the  effect  of  property 
variations  and  geometric  variations  associated  with  manufacturing  and  fabrica- 
tion. The  studies  included  in  Subsections  3.2  and  3.3  involve  axisymmetric 
variations  and  use  the  axisymmetric  finite  element  code.  In  Section  3.4 
asymmetric  effects  are  investigated  by  use  of  the  two-dimensional  finite 
element  model  shown  in  Figure  1-5.  In  this  latter  investigation  the  effect 
of  asymmetric  variations  on  the  Imbalance  of  the  flywheel  are  considered. 

3.2  VARIATIONS  IN  FIBER  VOLUME 

Here  we  consider  the  effect  on  the  flywheel  stresses  and  efficiency  of 
axisymmetric  variations  in  the  fiber  volume  ratio  V^.  Two  specific  cases  are 
considered!  (1)  the  case  of  one  constant  value  of  Vf  for  the  radial  layers 
and  another  constant  value  for  the  hoop  layers,  and  (2)  a case  of  a local 
variation  of  Vf  within  the  hoop  layer  alone. 

For  both  of  these  cases  we  assume  the  following  properties  for  the 
Kevlar  49  fiber  and  epoxy  matrix: 

Modulus  of  fiber  • Ef  - 20  x 10®  psl 

Modulus  of  matrix  • E^  • .5  x 10^  psi 

Fiber  weight/volume  ■ Pfg  ■ .053  lb/in3 

Matrix  weight/volume  - pn,g  - .0426  Ib/in^ 

The  Kevlar  49  properties  are  those  given  by  duPont.*  Typical  matrix  proper- 
ties are  used.** 

The  mechanical  properties  for  a composite  are  formulated  by  using  the 
rule  of  mixtures,  namely: 

Pg  - Pfg  »£  +P„g 

The  resulting  values  of  E as  a function  of  Vf  are  given  in  Figure  3-1.  It 
should  be  pointed  out  that  for  Vf  > .59  one  has  the  reference  properties 
El  « 12  X 10^  psl  and  pg  • .0487.  These  correspond  to  the  properties  cited 
in  Section  1. 

Case  1:  Different  Vf's  for  Hoops  and  Radisls 


In  this  case  we  fix  Vf  ••  .59  for  the  radial  reinforcements  and  vary  Vf 
for  the  hoop  reinforcements  over  the  range  Vf  ■ .45  to  Vf  ■ .70.  Thus  we 

•Hunter,  R.  L.,  "Cherectaristics  and  Uws  of  Kevlar  49  Hi(^  Modulus  Organic  Fiber,"  duPont  Report,  revised  9/18/73. 
••Ashton,  J.  E.,  Halpin,  J.  C.,  and  Petit,  P.  H.,  "Primer  on  Compotite  Materials:  Analysis"  Technomic  Publishers,  1969. 
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are  assigning  the  radials  the  reference  properties  and  we  are  varying  the 
properties  of  the  hoops.  Results  were  obtained  by  parametrically  running 
the  axisymmetric  finite  element  wheel  code.  Stress  distributions  obtained  in 
this  case  are  shown  in  Figure  3-2. 

In  order  to  determine  the  effect  of  Vf  on  the  wheel  efficiency  it  is 
necessary  to  postulate  the  dependence  of  the  strength  on  Vf.  This  was  de- 
termined on  the  basis  of  an  allowable  strain  criterion.  Consequently  the 
0^^  strength  is  assumed  to  be  proportional  to  the  modulus,  which  is  shown 
in  Figure  3-1. 

In  order  to  express  these  parametric  results  in  terms  of  energy  density 
or  a shape  factor  Kg  a problem  in  definitions  arises.  If  we  solve  the 
equation  in  Section  1.2  for  Kg  we  get: 

Kinetic  Energy  pg 
® Weight  a 

The  problem  is  to  decide  what  values  of  pg  and  a to  use  for  Case  1 in  which 
two  different  sets  of  properties  occur.  The  choice  used  in  this  work  is  to 
takepg  and  a as  the  reference  values,  i.e.,  those  for  Vf  > .59.  Having  done 
this,  all  variations  of  kinetic  energy  per  unit  weight  with  Vf  are  reflected 
in  corresponding  variations  of  Kg.  The  results  calculated  on  this  basis  are 
shown  in  Figure  3-1  which  shows  Kg  decreasing  as  one  moves  away  from  Vf  • .59, 
the  volume  fraction  for  which  the  wheel  design  is  optimized.  As  might  be  ex- 
pected, the  figure  shows  that  Kg  drops  more  rapidly  as  Vf  is  decreased.  The 
results  shovm  in  Figure  3-3  for  and  Oj.^  as  functions  of  Vf  do  not  indicate 
very  great  variations  within  the  range  Vf  ■ .45  to  .70,  and  are  not  expected 
to  significantly  limit  the  design. 

Case  2 1 Local  Variations  of  Vf  in  Hoops 

In  this  case  we  consider  the  effect  of  a local  axisymmetric  variation  in 
the  hoop  layer.  The  postulated  variation  of  Vf  is  shown  in  Figures  3-4  and 
3-5  together  with  the  stress  from  the  axisymmetric  analysis.  Also  shown  in 
these  plots  are  the  relative  changes  in  strength  postulated  on  the  basis  of 
strain  to  failure.  The  resulting  shape  factor  Kg  la  .447,  9 percent  down 
from  .474  of  the  original  design. 

3.3  VARIATIONS  IN  HOOP  AND  RADUL  GEOMETRIES 

Three  different  types  of  variations  in  the  geometry  of  the  hoop  and 
radial  reinforcements  were  studied: 

1.  Variations  of  hoop  thicknesses, 

2.  A .010  inch  change  in  radial  width, 

3.  A +.010  inch,  -.010  inch  combination  change  In  radial  width. 
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tionai  stress 


Variations  of  Hoop  Thickness 


Results  obtained  by  varying  the  hoop  thickness  are  shown  In  Figures  3-6 
through  3-8.  Figure  3-6  shows  the  changes  In  stress  distribution  that  occur 
as  the  hoop  thicknesses  are  perturbated  away  from  the  design  configuration. 

The  figure  Indicates  that  as  the  hoops  are  thickened  (+20%)  the  radlals  become 
overstressed  and  the  hoops  understressed.  When  the  hoop  thickness  Is  reduced 
(-20%)  the  reverse  occurs:  the  hoops  are  overstressed  and  the  radlals  become 
understressed.  These  results  clearly  show  the  optimal  feature  of  the  un- 
perturbed (0%)  original  design  thickness  combination.  The  corresponding 
changes  In  shape  factor  Kg  are  shown  In  Figure  3-7.  A 20  percent  reduction 
In  hoop  thickness  results  In  a 19  percent  decrease  In  Kg;  a 20  percent  In- 
crease In  hoop  thickness  lowers  Kg  by  only  10  percent.  The  wheel  appears  to 
be  more  sensitive  to  decreases  In  hoop  thickness  because  of  the  tendency  for 
Mgh  hoop  stresses  to  occur  at  the  Inside  edge  of  the  wheel. 

It  should  be  pointed  out  that  these  calculations  apply  equally  well  to 
changes  In  radial  thickness,  that  Is,  the  results  are  dependent  on  the  ratio 
of  hoop  to  radial  thickness. 

In  Figure  3-8,  the  corresponding  0^2  radial  Interlaminar  shear 

stresses  are  plotted.  The  special  nondlmenslonal  form  of  stress  shown  In  the 
figure  relates  the  shear  stress  to  the  maximum  stress  so  that  the  significance 
of  the  shear  stress  changes  can  be  evaluated.  We  presume  that  the  wheel  will 
be  spun  up  to  a speed  (1  for  ;4ilch  the  maximum  or  ^0  stresses  reach  the 
material  allowable  — consequently  the  ratio  <'rz/<^max  Is  a key  parameter  and 
Is  used  In  representing  the  results  in  Figures  3-8  and  3-11.  The  ordinates 
of  these  figures  are  expressed  In  terms  of  a shear  stress  parameter  (v^z’^)/ 
‘'max  *hf).  Figure  3-8  shows  that  the  effect  of  changing  hoop  thicknesses  of 
±20%  does  not  substantially  Influence  the  ‘'fz  shear  stress  parameter. 

.OlO-lnch  Change  In  Radial  Width 

In  the  reference  design  the  radial  reinforcements  are  .048  Inch  wide  on 
the  inside  (r  ■ 1.625"),  .191  Inch  wide  at  r » 4.85  Inch  s .a  .045  inch  wide 
at  the  outside  r " R - 9.75  inch.  We  now  consider  two  <. - yP-j  one  with  .010 
inch  added  to  the  radial  width,  and  one  with  .010  Inch  a<' ' ''..cted  from  the 
radial  width.  The  resulting  stress  distributions  and  : .ting  change  In 

shape  factors  are  shown  In  Figures  3-9  and  3-10.  The  distributions  are  very 
similar  and  there  Is  a small  change  In  Kg  of  approximately  8 percent  for  the 
+.010  and  -.010  changes  In  radial  width.  The  dominant  effect,  however,  is  a 
change  in  maximum  a^z  Interlaminar  shear  stresses  because  from  a percentage 
point  of  view  the  .010  inch  change  has  the  greatest  effect  near  the  Inside 
and  outside  edges  of  the  wheel  where  the  shear  stresses  are  greatest.  The 
changes  In  shear  stress  parameters  are  shown  In  Figure  3-11,  which  shows  an 
Increase  of  14  percent  for  the  radial  which  Is  .010  Inch  narrower. 

+.010  Inch.  +.010  Inch  Combined  Change  In  Radial  Width 

A combined  case  was  run  In  which  the  width  of  the  outer  half  (larger  r) 
of  the  radial  was  reduced  by  .010  Inch  and  correspondingly  Increased  on  the 


-42- 


Figure 


Local  radiut/outside  radius,  r/R 


i.  Sensitivity  of  stresiet  to  hoop  layer  thickness 
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Figure  3-9.  Stress  distributions  vtrsus  redial  width  variations  for  +.010  and  -.010  inch  changn 
in  radial  width. 
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Figure  3-10.  Smsitivity  of  shapo  factor  to  radial  width. 
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Inside  half  (smaller  r).  The  result  was  Kg  » .429,  a 9 percent  reduction 
from  the  design  case.  The  shear  stress  parameter  (ffrz*R)/(ffinax’hr)  14  per- 
cent higher  In  this  case. 

It  should  be  pointed  out  that,  for  the  +.010  Inch,  -.010  Inch  and  com- 
bined ±.010  Inch  cases,  that  the  maximum  shears  are  not  greatly  affected. 

3.4  EVALUATION  OF  ASYMMETRIC  PROPERTY  VARIATIONS 

Asymmetric  property  variations  were  evaluated  using  the  finite  element 
grid  shown  In  Figure  1-5,  The  model  "smears  out"  the  radlals  In  the  circum- 
ferential direction.  Two  cases  of  as3m)metrlc  variations  of  the  modulus,  as 
shown  In  Figure  3-12,  were  analyzed.  In  Case  A the  modulus  of  the  radlals  Is 
varied;  In  Case  B the  modulus  of  the  hoops.  The  change  In  stresses  for  these 
cases  Is  as  follows: 


Case 

Change  In  modulus 

Change  In  maximum  stress 

A 

±5% 

+3.9% 

B 

±1% 

+.1% 

A major  reason  for  carrying  out  these  asymmetric  property  variation 
studies  Is  to  compute  the  Imbalance  due  to  deformation  of  the  wheel.  To  do 
this  we  evaluate  the  shift  In  centroid  of  the  wheel  (relative  to  Point  A of 
Figure  3-12)  that  would  exist  when  Omax  “ 200  ksl.  A typical  wheel  weight  of 
16.86  pounds  was  used  In  the  calculations.  The  results  presented  In  tabular 


form  appear 

below: 

Case 

Change  In  modulus 

Distance  of 
centroid  shift 
(Inch) 

Imbalance  for  16.86 
pound  wheel 
(Inch-grams) 

A 

±5% 

.00416 

36.19 

B 

±i% 

.000(04 

6.9 

The  resulting  Imbalances  shown  above  are  considerably  In  excess  of  the  Im- 
balances that  remain  after  balancing  the  wheel  which  are  of  the  order  of 
1 Inch-gram.  It  should  be  recognized  that  the  ±5%  and  ±1%  variations  used 
In  the  study  are  highly  unlikely  to  occur  because  they  represent  averages 
of  all  the  deviations  In  modulus  through  the  thickness  as  well  as  In  the 
In-plane  directions.  The  variations  used  In  this  study  are  conservative 
estimates  of  asymmetries  In  an  Individual  layer.  Random  layer  orientation 
during  assembly  and  Inter-layer  load  sharing  during  operation  will  tend  to 
distribute  displacement  more  evenly.  The  method  of  fabricating  the  hoops  Is 
more  likely  to  generate  axlsymmetrlc  rather  than  asymmetric  variations  In 
modulus.  Similarly,  the  radlals  are  stamped  out  of  sheets  of  Kevlar  tape 
and  are  expected  to  be  very  regular. 
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KEY 

Eg  ■■  Modulus  of  hoops 

Ef  Modulus  of  radials 

E « Referanca  modulus,  12.0  X 10^  psi 


Figurt  3-12.  2-0  Modalt  us«d  to  •valuato  atymmetric  proptrty  variations. 
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SECTION  4 


ANALYSIS  OF  STATIC  SIMULATION  TEST 


4.1  INTRODUCTION 

The  purpose  of  this  investigation  Is  to  attempt  to  design  a static  test 
that  will  simulate  the  flywheel  stress  conditions  that  would  be  encountered 
under  inertial  loading  In  practice.  Such  a test  would  of  course  greatly  ex- 
pedite the  development  of  bi-directionally  reinforced  fl3rwheels  by  avoiding 
the  problems  associated  with  dynamic  testing  and  by  allowing  test  specimens 
to  be  thinner  and/or  smaller  than  a full  scale  flywheel  and  thus  save  on 
fabrication  costs. 

As  a starting  point  we  will  seek  to  develop  a test  configuration  that 
will  simulate  the  stress  components  ‘'rd  shown  In  Table  l-l. 

The  2-D  and  equivalent  axlsymmetrlc  results  shown  In  Table  1-1  are  very 
similar,  and  we  shall  pick  as  a starting  point  those  given  by  the  equivalent 
axlsymmetrlc  analysis  because  It  was  used  originally  to  design  the  wheel. 
Expressing  these  stresses  In  relative  terms  gives  the  following: 


O^loQ  ■ 1.0 

(4-1) 

a loa  ■ .00777 
rz  ^ 

(4-2) 

OQ  » .01903 

(4-3) 

It  should  be  pointed  out  that  these  ratios  do  not  occur  simultaneously  In  the 
original  flywheel  because  a-^  and  oq  do  vary  with  the  local  radius  r.  However, 
as  a first  cut  at  the  simulation  problem  It  seems  reasonable  to  seek  a simula- 
tion as  close  as  possible  to  the  relative  values  cited  above. 

The  strategy  used  to  find  a configuration  which  simulates  the  three  stress 
ratios  Is  similar  to  that  used  to  design  the  original  flywheel.  First  equiva- 
lent "smeared  out"  axlsymmetrlc  analyses  and  other  simplified  analyses  are 
carried  out  to  obtain  geometric  configurations  approximately  simulating  the 
desired  stress  ratios.  As  In  the  flywheel  design,  the  E90  ■ 0 assumption  Is 
carried  along  to  simulate  a crazed  condition.  After  this,  a two  dimensional 
finite  element  analysis  Is  carried  out  to  verify  the  design. 

4.2  GEOMETRIC  PARAMETERS 

The  geometric  parameters  available  In  the  simulation  test  shown  in  Fig- 
ure 4-1  are: 

N^  The  number  of  radlals  per  layer 

Rj  , R ^ The  inside  and  outside  radii  of  the  hoops 

in’  out 
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I 


The  radial  and  hoop  thicknesses  (actual) 

y (r)  The  portion  of  the  circumference  at  a given  radius  r 

occupied  by  radials. 

The  ratio  of  internally  applied  load  to  externally  applied  load  is  also  a 
parameter.  However,  for  the  Initial  part  of  this  Investigation  the  Internal 
pressure  Is  assumed  to  be  zero. 


For  analytical  purposes  It  Is  convenient  to  Introduce  two  more  parameters 
as  follows: 

h^  (r)  "smeared  out"  radial  thickness,  averaged  with  respect  to 

variations  with  $ 

hg  (r)  "smeared  out"  hoop  thickness,  analogous  to  h^  (r) . 

In  this  work,  we  consider  only  the  case  of  uniformly  thick  hoop  layers  such 
that  ■>  H^. 

A. 3 SIMULATION  OF  - 1 

Ideal  Shape 

For  the  equivalent  "smeared  out"  axisymmetric  model  with  E90  *0,  it  Is 
possible  to  solve  for  fi^Cr),  the  equivalent  radial  thickness  for  which  o^^/oq  ■ 
1 throughout  the  region  < r < Rout* 

Under  these  assumptions  the  equation  of  equilibrium  for  the  case  ^ o q 
and  n • 0 (no  Inertial  force)  Is; 


-or  + ff  h,+  ar  -7-  = 0 
dr 


(A-A) 


where  Is  the  thickness  of  the  hoop  layer,  which  Is  constant  with  respect 
to  r. 


The  solution  of  this  equation  subject  to  a force  free  inner  boundary 
condition  of  a*hj  ■ 0 which  becomes  hr  ■ 0,  is  the  following: 


f 


A plot  of  this  equation  Is  given  In  Figure  A-1. 

Unfortunately,  this  Ideal  shape  leads  to  Infinitely  high  interlaminar 
shear  stresses,  o between  the  hoops  and  radials  at  the  Inside  edge  (r  « 
Rin) . This  comes  about  in  the  following  way.  First,  from  Equation  (A-5)  one 
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f 


I?'- 

has  hr  - 0 when  r « Rin.  Second,  from  the  definition  of  the  "smearing  out", 
the  width  of  a radial  (i.e,,  in  0 direction),  W(r) , is  given  by: 


W(r)  = 


\ 2 


TTt 


NrHr 


hr(r) 


Consequently,  at  the  inside  edge  (r  » Rin)»  the  radial  width  is  zero.  Thus, 
as  one  approaches  the  inside  edge  (r  -» R^n)  the  area  for  load  transfer  between 
hoop  and  radial  reinforcements  approaches  zero  and  infinitely  high  interlaminar 
shear  stresses  arise.  This  effect  can  also  be  deduced  from  Equation  (4-14). 

Solution  for  Constant  h^ 

7 

This  problem  was  solved  in  the  case  of  the  wheel  by  increasing  hr  near 
the  inside  edge,  i.e.,  by  adding  "tabs."  In  the  simulation  we  anticipate  the 
case  of  a relatively  narrow  ring  (RoufUln  Is  small)  and  therefore  seek  a 
configuration  like  that  shown  in  Figure  4-2  where  the  radials  occupy  a constant 
angular  portion  of  the  wheel.  This  analytical  model  can  be  solved  by  solving 
the  equation 


do 


f * 


where  now  hr  ^s  a constant  and  (^r  and  oq  are  given  by 


(4-6) 


w 


(4-7) 

(4-8) 


The  solution  obtained  by  substituting  oq  and  Oj.  into  the  equilibrium  equation 
and  invoking  the  boundary  condition  ar."  0 at  r ► Rin  is  as  follows: 


[f- 

L\''« 

lUt/  V^OUt  / 

T] 

a a 

[(£)■■■ -1 

(4-9) 
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where  o is  the  value  of  at  the  outside  radius  r » Rout  given  by 


(4-11) 

It  is  desired  to  pick  a ratio  hj./hg  such  that  the  maximum  values  of 
and  oq  are  equal.  A plot  giving  the  desired  ratios  h^/h^  as  determined  from 
the  above  formulas  is  given  in  Figure  4-3. 

Hoop  Bending  Effects 

If  the  radlals  shown  in  Figure  4-2  are  too  far  apart,  then  the  hoop  ring 
will  undergo  bending  and  the  values  of  and  oq  will  deviate  from  those  pre- 
dicted from  the  equivalent  "smeared  out"  axlsymmetrlc  analysis.  Although  this 
problem  Is  negligible  In  the  actual  flywheel  where  (nuinber  of  radlals)  and 
Rout~Rln  large  It  Is  of  primary  concern  In  a simulation  test  where  It  Is 
desirable  to  keep  as  low  as  possible  for  practical  and  cost  reasons. 

This  effect  was  evaluated  parametrically  by  solving  Flugge's  ring 
equations*  for  the  case  of  a simplified  ring  model  where  the  radial  loads 
are  assumed  to  be  applied  as  point  loads  at  the  position  of  each  radial 
centerline.  This  analysis  Is  expected  to  overestimate  the  effects  caused 
by  bending  and  provide  a conservative  method  for  excluding  simulation  tests 
designs  with  Insufficient  numbers  of  radlals.  Applying  Flugge's  equations 
and  solving  for  the  moment  Mo  In  the  ring  at  the  point  of  load  application 
gives 


Mg  = — iQg  R (coo.  d - 1)  + Ng  R (6  - 0 sin  0) 

« L 

where 

D >■  polar  angle  between  centerline  of  radial  and  middle  point  between 
radlals . 

R*  average  radius  of  the  hoop  ring 


(4-12) 


Flugge.W.,  Stresses  in  Shells, Sprinper-Verlag,  Berlin,  1960. 


Qg  > shear  force  at  centerline  of  radial  (equals  on  half  of  total  radial 
load) 

N = hoop  membrane  force, 
o 


It  was  found  that  Mg  Is  greater  In  absolute  magnitude  than  the  bending  moment 
midway  between  radlals.  Computing  the  bending  stress  due  to  Mg  and  comparing 
It  to  the  hoop  stress  associated  with  Ng  leads. to  the  design  plot  shown  In 
Figure  A-4.  The  asymptotic  formula  for  small  0 la 


lending 

‘^hoop 


A A A 

2 R - 


(4-13) 


It  can  be  seen  from  Figure  4-4  that  If  a typical  case  of  Rout'^ln  “ ^ 
taken,  then  approximately  25  radlals  (Ny  ■ 25)  are  required  to  keep  the  bend 
stresses  below  5 percent  of  the  hoop  stresses. 


It  turns  out,  however,  that  the  restrictions  Implied  by  Figure  4-4  are 
less  severe  than  those  required  to  simulate  the  condition  o^q  Ioq  - .01903. 


4.4  SIMULATION  OF  - .00777 

The  derivation  used  to  compute  and  for  the  equivalent  "smeared 
out"  axlsymmetrlc  analysis  Is  shown  schematically  In  Figure  4-5.  The 
equations  are  obtained  by  taking  free  body  diagrams  of  EFCD  and  ABCD, 
respectively.  The  assumption  Ego  ■ 0 precludes  the  possibility  of  any 
radial  stresses  a,,  on  the  surfaces  FC  or  ED  of  the  hoop  layers.  These 
surfaces  can  contain  very  small  amounts  of  shear  stress,  but  as  Indicated 
by  the  more  detailed  2-D  finite  element  analyses  they  act  almost  perpen- 
dicular to  the  radial  direction  used  In  summing  forces  and  are  negligible. 

The  factor  of  2 in  the  equation  for  average  Is  caused  by  the  fact  that  the 
load  Is  transferred  from  two  sides  of  each  hoop  layer.  As  in  the  design  of  the 
flywheel,  this  has  serious  implications  for  layers  adjacent  to  the  outside 
edge  (Z  ■ flywheel.  In  the  flywheel  design,  an  extra  wide 

set  of  radlals  was  used  Just  Inside  the  outermost  hoop  layer  to  maintain 
similar  levels  of  cr^z  stress  near  the  outside.  For  the  simulation  it 
appears  more  appropriate  to  put  a set  of  radlals  one  half  the  usual  thickness 
on  the  outside  so  that  the  stresses  do  not  vary  from  layer  to  layer. 

Working  out  the  geometry  indicated  In  the  equations  of  Figure  4-5  leads 
to  the  following: 


rz 


ave. 


(4-14) 
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N » Number  of  radials/layer 


50 
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Figure  4-4.  Percent  chenge  in  hoop  streu  cauied  by  local  bending. 
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I 


pJJ^rdV 


Summing  forcM  in  local  radial  direction!  gives 

a„  x2xAREA(AEFB)-^o-DChp-Oa 
rtave  ‘ w w 

®rd  • ^’6  • AB  ■ (pO^  r^  - dp)  • hg  • DC 
For  static  simulation  R - 0 

97-731 

Figura  4-6.  Sch«matie  of  darivation  of  Ofj,  and  Of^  ttratsai . 
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where 


r Is  the  local 
y is  the  local 


radius 

fraction  of  the  circumference  occupied  by  radlals. 


4.5  SIMULATION  OF  o-cQI%  " .01903 

Similar  application  of  the  equations  shown  in  Figure  4-5  gives 


(4-15) 

f 

where  N is  the  number  of  radials. 
r 

As  we  shall  see  in  the  subsequent  section,  this  last  equation  places  the 
most  severe  restraint  on  the  minimum  acceptable  number  of  radlals,  N^. 

4.6  SELECTION  OF  A TYPICAL  SIMULATED  FLYWHEEL  CONFIGURATION 

We  next  apply  the  relationship  developed  in  Sections  4.1  through  4.5  to 
define  a simulation  test. 

In  order  to  satisfy  Equation  (4-15)  a tradeoff  must  be  made  between  the 
dual  evils  of  too  many  radlals  on  the  one  hand  and  too  large  a y,  the  relative 
area  occupied  by  radials,  on  the  other.  Our  mutual  selection,  satisfying 
Equations  (4-15)  and  (4-3),  is: 

N^  - 48 

y - .7092 

If  we  select  ■ 1.625  - same  as  flywheel,  then  the  circumferential  distance 
between  radlals  centerlines  is  .2127.  The  areas  between  radial  centerlines 
are  then  kept  roughly  square  by  selecting  - Rj^jj  - .25  inch.  For  this 
combination  Equation  (4-13)  indicates  an  acceptable  deviation  from  pure  hoop 
stress  of  6 percent.  For  this  combination  of  Rin/l^out'  4-3  Indicates 

hr/h^  ■ .125. 

We  next  use  Equations  (4-14)  and  (4-2)  to  find  h<?  » .0179.  Then  for  hj. 
we  have  .125  x .0179  ■ .0022  inch.  The  actual  radial  thickness  * 

.00310.  In  summary,  the  resulting  simulated  configuration  is; 


Geometry  Selected  for  Simulation  Test 


Nr 

■ 49  radials /layer 

y 

- ,7092 

hn 

- 1,625 

*^out 

- 1.875 
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.00310  actual  radial  thlckneas 

h - .0022 

r 

he  - .0179 

If  It  Is  found  impractical  to  fabricate  the  values  of  Hj-  and  cited 
above,  one  can  merely  scale  up  ^out*  ^r  until  a more  suitable 

combination  is  found,  without  altering  the  accuracy  of  the  simulation. 


-62- 


SECTION  5 


TESTING  FOR  DEGRADED  PROPERTY  EFFECTS 


5.1  MECHANICAL  EVALUATION 

A biaxial  tensile  test  was  devised  whereby  the  Kevlar  0°  - 90°  lay-up 
could  be  stressed  in  both  directions  simultaneously.  Subsequent  to  this 
biaxial  loading  the  stressed  composite  at  the  junction  or  overlap  region  was 
machined  into  samples  for  evaluation  of  shear  properties  after  loading. 

The  method  of  testing  in  this  series  is  a sandwich  beam  using  a honey- 
comb core  between  an  aluminum  skin  and  a composite  (test  piece)  skin.  The 
loading  method  Is  essentially  a double  beam  flexure  test  using  four  point 
loading  on  each  leg.  The  schematic  of  this  test  is  given  in  Figure  5-1.  The 
samples  are  loaded  at  particular  stress  levels  to  simulate  various  flywheel 
speeds,  l.e.,  50,  60  and  IQX  of  ultimate,  and  then  off-loaded.  The  composite 
material  is  then  removed  from  the  sandwich  structure  and  samples  are  fabri- 
cated from  the  areas  at  the  0°/90°  overlaps.  The  samples  are  designed  to 
obtain  shear  properties  after  biaxial  loading.  The  tests  selected  are  45° 
tension  for  determining  Young's  modulus  and  Poisson's  ratio  (which  could  then 
be  converted  to  shear  modulus)  and  tensile  and  shear  stress;  and  horizontal 
(short  beam  shear)  to  determine  interlaminar  shear  stress.  Figure  5-2  shows 
the  areas  and  orientations  of  the  test  samples  removed. 

The  matrix  of  tests  to  be  performed  In  this  series  Is  given  In  Table  5-1. 


TABLE  6-1.  0°/90“  BEAM  TEST  (BIAXIAL  TENSION) 


Span 

Load  point  separation  Load  level 

Number  of  replicates 

(Inch) 

(inch) 

(Z) 

20 

4 

100 

2 

20 

4 

80 

2 

20 

4 

70 

2 

20 

4 

60 

2 

SUBSEQUENT  COUPON  TESTS  ON  COMPOSITE 

OVERLAP 

No.  of  replicates 

Test 

(each  beam) 

Data  generated 

45°  tensile 

1 

45°  , G.  E, 

45° 

2 

(Interlaminar) 
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5.2  COMPOSITE  LAY-UP 


Since  there  was  an  abundance  of  9 Inch  long  uni-directional  prepreg 
Kevlar  epoxy  on  hand  at  the  outset  of  this  test  series,  and  since  there  was 
a waiting  period  for  receipt  of  any  roll  goods  which  were  on  order,  it  was 
decided  to  go  along  and  fabricate  the  first  two  beams  with  9 inch  strips 
layed  up  with  overlap  and  butt  joints.  The  remaining  beams  in  the  series, 
however,  are  made  with  22  inch  strips  of  uni-directional  Kevlar  layed  up  at 
0°  and  90°  alternating  angles.  Also,  the  first  two  beams  contain  2 plies  of 
composite  in  each  direction  for  a total  of  4 plies  whereas  the  remaining 
samples  are  constructed  with  3 plies  in  each  direction.  Figure  5-3  is  an 
attempt  to  depict  the  butt  and  overlap  joint  sequence  herewith  described. 

5.3  REDUCTION  OF  DATA 

At  the  time  of  this  writing,  one*  test  has  been  carried  out  at  52  per- 
cent of  ultimate  load.  From  these  tests  the  following  average  properties 
were  obtained: 

E45O  - .99  X 10^  pal 

^ultimate  " ^750  psl 

The  next  step  is  to  derive  from  these  quantities  the  quantities  Gj^2  ‘'12 
allowable.  The  modulus  Gyi  obtained  from  the  transformation  laws  and  is 


which  gives  Gj^2  * *25  x 10^.  The  state  of  stress  at  a 45°  orientation  from 
the  direction  of  loading  is  given  by 

a 

<7j  = <72  = «'12  = y 


Tsai-Hill  Failure  criterion  then  gives 


= 1 


*A  second  biaxial  test  failed  and  prevented  measurement  of  the  degraded  modulus. 
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97733 


Figurt  5-3.  Butt  jointt  and  overlap  schama  for  8 -inch  long  strips. 
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which  reduces  to: 


Since  for  Kevlar  49  the  undegraded  shear  strength  S is  much  lower  than  the 
0°  strength  Y we  have 

S = = 5750 


Tabulating  the  results  gives 


Quantity 

Undegraded 

Degraded 

Ratio  of  degraded/undegraded 

Gi2 

.3  X 10^  psi 

.25  X 10^ 

72. 8X 

0^2  allowable 

7.9  ksl 

5.75  ksi 

82. 5X 

These  results  are  then  Introduced  into  the  analysis  of  Section  2 to  evaluate 
the  effects  of  degraded  properties  on  flywheel  performance. 
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SECTION  6 


CONCLUSIONS 


Structural  analyses  and  experiments  have  been  carried  out  to  further  the 
development  of  bi-directionally  reinforced  flywheels.  A brief  summary  of  the 
findings  are  presented. 

1.  The  analysis  Indicates  that  for  the  biaxial  stress  levels  used  in 
the  testing  that  shear  modulus  and  strength  degradations  would  not 
give  rise  to  negative  stress  margins.  In  addition,  the  analytical 
results  show  that  reductions  in  in-plane  shear  modulus  of  up  to  90 
percent  have  a small  effect  on  the  maximum  0°  stresses. 

2.  The  efficiency  of  a bi-directionally  reinforced  flywheel  is  sensitive 
to  variations  in  fiber  volume,  and  variations  in  hoop  and  radial  re- 
inforcement geometry.  The  relationship  between  variations  in  effi- 
ciency and  variations  in  property  and  geometric  inputs  differ  from 
case  to  case  and  are  presented  in  detail  in  the  report.  Asymmetric 
variations  in  material  properties  can  also  lead  to  substantial  wheel 
Imbalances  which  depend  on  the  degree  on  nonuniformity  present. 

3.  Static  simulation  tests  appear  feasible  but  must  include  a sufficient 
number  of  radlals  to  prevent  premature  shear  failures. 

4.  Blaxially  degraded  property  effects  have  been  obtained  from  crossbeam 
sandwich  tests  and  subsequent  45°  tensile  tests  and  the  resulting 
data  input  into  the  analysis. 
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